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Hyperdoping has emerged as a promising method for designing semiconductors with unique optical 
and electronic properties, although such properties currently lack a clear microscopic explanation. 
Combining computational and experimental evidence, we probe the origin of sub-band gap optical 
absorption and metallicity in Se-hyperdoped Si. We show that sub-band gap absorption arises from 
direct defect-to-conduction band transitions rather than free carrier absorption. Density functional 
theory predicts the Se-induced insulator-to-metal transition arises from merging of defect and 
conduction bands, at a concentration in excellent agreement with experiment. Quantum Monte 
Carlo calculations confirm the critical concentration, demonstrate that correlation is important to 
describing the transition accurately, and suggest that it is a classic impurity-driven Mott transition. 



Of all the experimentally measurable physical proper- 
ties of materials, the electronic conductivity exhibits the 
largest variation, spanning a factor of 10^^ from the best 
metals to the strongest insulators^. Over the last century, 
the puzzle of why some materials are conductors and oth- 
ers insulators, and the mechanisms underlying the trans- 
formation from one to the other, have been carefully scru- 
tinized; yet even after such a vast body of research over 
such a long period, the subject remains the object of 
controversy. In 1956, Mott introduced a model for the 
insulator-to-metal transition (IMT) in doped semicon- 
ductors, in which long-ranged electron correlations are 
the driving forc^. Hyperdoping (doping beyond the sol- 
ubility limit) creates a new materials playground to ex- 
plore defect-mediated IMTs in semiconductors. In this 
letter, we identify a defect-induced IMT in silicon hyper- 
doped with selenium concentrations exceeding 10^^ cm~^ 
(compared to the equilibrium solubility limit^ of about 
10-^^ cm~^) and we explore the detailed nature of the 
transition with both experiment and computation. We 
find that the IMT is largely driven by electron correla- 
tion and most resembles a classic impurity-driven Mott 
transition. Additionally, we find that the high density 
of Se present at the IMT yields direct optical transitions 
and an absorption coefficient in excellent agreement with 
the measured sub-band gap optical properties^. 

Hyperdoping is currently being used to engineer new 
materials with unique and exotic properties. Silicon 
hyperdoped with chalcogens exhibits strong sub-band 
gap absorption down to photon energies as low as 0.5 
eVMHI^ sparking substantial recent interest in applica- 
tions such as infrared detection and intermediate band 
photovoltaics^ The successful fabrication of rectify- 
ing junctionals) photodiode^^^^^^ using S and Se hy- 
perdoped silicon appears to justify such interest. While 
isolated S and Se dopants are well-established deep dou- 
ble donors in silicoiP^, the enhanced optical properties 
of hyperdoped silicon (in which these chalcogenic impu- 



rities are present at much higher concentrations) are not 
yet well understood. Further, unlike the prototypical sys- 
tem of phosphorus-doped silicon for whic h the IMT has 
been extensively studied and characterizecP^^, there are 
very few studies of an IMT resulting from deep defects 
such as chalcos ensP. 

We prepared Se-doped silicon (Se:Si) samples using ion 
implantation followed by nanosecond pulsed-laser melt- 
ing (PLM) and rapid resolidification. The PLM pro- 
cess enables chalcogen doping with concentrations ex- 
ceeding 1% atomic; such samples exhibit unexplained op- 
tical properties including broad, featureless absorption of 
photons with energy lower than the band gap of silicorP. 
Silicon substrates (boron doped, p ~ 25 l^-cm) were ion 
implanted with Se to nominal doses of 3(10^^) and 10^^ 
cm^ using an ion beam energy of 176 keV. The implanted 
samples were exposed to four laser pulses (fiuences of 1.7, 
1.7, 1.7 and 1.8 J cm~^). This fiuence regimen results in a 
slightly shallower dopant profile, and higher peak Se con- 
centration, than reported previousl}/^^. The Se-rich layer 
is crystalline, extends approximately 350 nm from the 
surface, and is electrically isolated from the p-type sub- 
strate by the rectifying junction formed between the two. 
The Se concentration-depth pr ofile w as measured via sec- 
ondary ion mass spectrometrjEEHl. Sample preparation 
and measurement proceeded as described previousljEl, 
with the Se:Si layer isolated using cloverleaf-mesa struc- 
tures. Conductivity was calculated from sheet conduc- 
tivity using the effective implantation depth d^ff^^. 

Since charge carriers in the insulating phase require 
thermal activation, rigorous experimental proof of an 
IMT lies in the measurement of nonzero conductivity as 
the temperature T approaches K. In Fig.[T] we plot the 
temperature-dependent conductivity over a temperature 
range 1.8-300 K for two samples exhibiting peak Se con- 
centrations of 1.4(10^^) and 4.9(10^°) cm"^. Despite the 
relatively small (~3.5x) difference in peak dopant con- 
centration, the conductivity varies by almost 5 orders of 
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FIG. 1. Temperature— dependent conductivity of Se— 
hyperdoped silicon. At peak Se concentrations of Sei:Si349, 
hyperdoped Si exhibits strongly temperature-dependent con- 
ductivity, indicative of the insulating phase. At higher con- 
centrations (SeiiSigg), the conductivity is comparatively in- 
sensitive to temperature down to T = 1.8 K, indicating a 
dopant-induced transition to the metallic state. Inset: Anal- 
ysis of density functional theory (DFT) calculations (shown in 
Fig. [2| indicate that the conduction and defect bands cross 
as the Se-dopant concentration increases. Conduction and 
defect band edges are demarked by DFT's Kohn-Sham eigen- 
values for the insulating systems, with energies referenced to 
the valence band maximum ^vbm- 



magnitude at T = 1.8 K. The dramatic difference in con- 
ductivity between two samples doped to similar levels, 
and the lack of significant temperature sensitivity as T 
approaches K in the more highly-doped sample, demon- 
strate a transition to metallic conduction at a concentra- 
tion between that of the samples presented here. For 
ease of subsequent comparison to computational results, 
we identify this concentration by the peak value of the 
of Se to Si ratio in these samples, thus between Sei:Si349 
and SeiiSigg. 

While the experimental evidence for an IMT is clear, 
it cannot indicate the origin of the transition. The mea- 
sured critical Se concentration is accessible, though, to 
both density functional theory (DFT) and the more ac- 
curate quantum Monte Carlo (QMC) methods. For shal- 
low donors in silicon, in contrast, the transition occurs 
at roughly 1 dopant per 10000 atomM The DFT^^ 
results presented here employ the PBE approximation^^ 
to the exchange correlation functional as implemented 
within the SIESTA packag^. The inner core electrons 
are represented by Troullier-Martins pseudopotentials, 
and the Kohn-Sham orbitals are represented via a lin- 
ear combination of numerical pseudo-atomic orbitals ex- 
panded in a triple-zeta with polarization gaussian basis 
set. The DFT predicted lattice constants and band gaps 
are, respectively, for silicon 5.48 A(exp: 5.43 A) and 0.70 
eV (exp: 1.1 eV), and for hexagonal selenium a=4.40 



A, c/a=1.16 (exp: a=4.37 A; c/a=1.14) and 0.85 eV 
(exp: 1.8 eV). For each doping concentration, we substi- 
tuted one silicon atom with a Se impurity and relaxed the 
atomic positions (so that atomic forces are <0.01 eV/A) 
and supercell lattice vectors (so that all stress tensor com- 
ponents are <2 kbar). We chose the substitutional Se 
configuration, as our electronic structure calculation^^ 
as well as others^^ indicate the substitutional site to be 
the minimum energy defect configuration. The QMC re- 
sults presented here are computed via fixed node diffu- 
sion Monte Carlo conducted with the QWalk cod^^, with 
trial wave functions constructed with a Slater-Jastrow 
form using SIESTA's DFT orbitals, variance-minimized 
Jastrow coefficients, and a time step of 0.01 au. The 
QMC energies are computed by averaging over twisted 
boundary conditions for all supercells. Defect formation 
energies for both DFT and QMC were computed using 
AEf = {Esei-.Sir, + I^Si) - {Esir,+i + M5e), wherein each 
atom's chemical potential ja is determined by the quasi- 
equilibrium of the doped silicon with SiSe2 chains, which 
may be present as early stage precipitates^. 

Using DFT, we first explore the electronic band struc- 
ture, the electronic density of states (DOS), and the de- 
fect formation energy for supercells of size n x n x n, 
for n = 2,3, ...,8. Using the 2-atom face-centered cu- 
bic unit cell, this sampling corresponds to systems of 
Sei:SeAr-i, for N = 1024,686,432,250,128,54, and 16, 
and defect spacing that increases uniformly from one su- 
percell to the next. The left-hand column of Fig. [2] shows 
the computed DOS for all concentrations (shaded por- 
tions denote occupied states). At the lowest concentra- 
tion (Sei:Siio23), an isolated defect peak of narrow band- 
width (~0.06 eV) is offset from the conduction band edge 
by ~0.12 eV (estimated from the Kohn-Sham eigenval- 
ues). This result is consistent with previous work^^. The 
defect band is completely occupied by the two extra elec- 
trons introduced by the Se impurity. Because the filled 
defect band is offset from the conduction band, there are 
no nearby empty states and these extra electrons cannot 
contribute to metallic conduction. Thus, the system is 
in the insulating state. The middle column of Fig. [2] il- 
lustrates the charge density of the defect level, plotted at 
the F-point and on the Si (111) plane. At the concentra- 
tion of Sei:Siio23 the defect level charge density is highly 
localized around the impurity. 

As the impurity concentration increases from 
Sei:Siio23 to Sei:Si249, the strongest effect evident in 
the DOS in Fig. |2] is the increasine; width of the defect 
peak, indicating increasing dispersion and interactions 
between neighboring defects (also apparent in the defect 
state charge densities in Fig. 2]). As defects become 
more closely spaced, stronger defect-defect interactions 
result in a decrease of the offset between the defect 
peak and the conduction band minimum in the DOS. 
At the concentration of Sei:Si249 the defect peak just 
touches the conduction band edge, and at Sei:Sii27 
the defect and conduction bands have merged. The 
disappearance of the offset between the defect peak and 
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FIG. 2. Evolution of the insulator— to— metal transition in selenium— hyperdoped silicon. Left column: the electronic 
density of states (referenced to the Fermi energy Ef) for a range of Se concentrations, illustrating that the defect- and 
conduction-bands merge at Sei:Si249. The y-axis is uniform for all concentrations, and the DOS is normalized so that each 
plot contains the same number of atoms; shaded regions indicate occupied states at T = K. Middle column: charge density of 
the defect state at the F-point, plotted on the Si (111) plane. Increasing dopant concentration increases interactions between 
neighboring defects, resulting in the eventual delocalization of the defect state. Right column: the calculated absorption 
coefficient arising from both (red) direct optical transitions (for pure and metallic Se:Si) and (blue) free carrier absorption (for 
metallic Se:Si), in comparison to pure silicon (black). The green line for Sei:Sii27 indicates to the experimental data of Ref. [4]. 



the conduction band edge at Sei:Si249 indicates that 
the dopant electrons are no longer bound to impurity 
sites. Instead, many low-lying conduction-band like 
states are available for charge transport without thermal 
activation, signifying the onset of the IMT. This band- 
crossing induced phase transition is also illustrated in 
Fig. [l] From Figs, [l] (inset) and [2) DFT predicts that 
the transition occurs at a concentration '^Sei: 81249; this 
density is in excellent agreement with experiment. 

For defect concentrations of Sei:Sii27 and higher, the 
defect band further decomposes as it merges with the 
conduction band, as shown in Fig. [2] and Fig. [l] Large 
amounts of charge between impurities is observed in the 



charge densities of Fig. |2j indicating significant delocal- 
ization of the defect states. The DFT analysis indicat- 
ing hybridization of defect states with the conduction 
bands is consistent with experimental results, which in- 
dicate that the charge carriers exhibit conduction-band 
like character. For example, we calculated the Hall mo- 
bility of the metallic sample in Fig. [T] by measuring the 
Hall voltage as well as conductivity. The Hall mobility 
of the metallic sample is 21 ± 2 cm^ V"^ s"^ at T = 4.2 
K, comparable in magnitude to that of metallic silicon 
doped with shallow donors (for w hich t ransport is known 
to occur in the conduction band)P^^. This high mobil- 
ity value suggests that conduction in the metallic Se:Si 
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FIG. 3. Defect formation energy vs. defect concentra- 
tion. Formation energy of a substitutional Se dopant in Si, 
computed via DFT and QMC. According to DFT, a peak ap- 
pears in the formation energy near the concentration Sei :Si249 
corresponding to the IMT. This peak is enhanced and shifted 
towards higher concentrations (Sei:Sii27) in the QMC calcu- 
lations, as is the phase transition. Also, the QMC correction 
to DFT is largest near the transition, suggesting that electron 
correlation effects are prevalent at these concentrations. 



sample arises from delocalization of Se electrons into the 
conduction band. Thus both experiment and calculation 
support the idea that the IMT occurs due to a band- 
crossing of the Se defect states with the conduction band. 

In the right-hand column of Fig. [2j we show contri- 
butions to the absorption coefficient arising from direct 
optical transitions calculated using our DFT results. In 
addition, contributions arising from free carrier absorp- 
tion are shown, obtained using the Drude description and 
the experimentally measured free carrier concentration 
for the metallic Se:Si sampl^^. Figure [2] demonstrates 
that the defect levels predicted by our DFT calculations, 
and the resulting direct optical transitions and absorp- 
tion coefficient, yield excellent agreement with both the 
magnitude and shape of the measured sub-band gap op- 
tical absorption in Se-hyperdoped silicoiP^. Free carrier 
absorption, alternatively, does not become comparably 
important except at photon energies lower than 200 meV. 

We now consider the underlying nature of the IMT 
itself. In 1956, Mott introduced a model for an IMT 
in doped semiconductors driven by long-ranged elec- 
tron correlation^, proposing that the delocalization of 
the electronic wave function occurs through long-ranged 
electron-electron interactions that screen the binding 
field of the impurity potentiaP^. In the extensively 
studied semiconductor systems (such as P:Si), the mea- 
sured critical concentration for metallic onset coincides 
well with concentration predicted by the Mott criterioiP. 
Alternative IMT mechanisms have been identifiecF', in- 
cluding transitions arising from static disorder (Anderson 
mechanism), local electron-electron correlations (Mott- 
Hubbard mechanism), and strong electron-phonon cou- 



pling; in real systems, more than one mechanism may 
contribute simultaneously. Definitive proof that the ob- 
served IMT is driven by long-ranged electron correla- 
tion is difficult, but QMC can establish the relative im- 
portance of electron correlations on the defect forma- 
tion energies in both insulating and metallic systems. 
DFT treats electron correlation in an approximate man- 
ner (correlation refers here to contributions to the total 
energy beyond the independent electron approximation, 
including the exchange contribution), and describes to- 
tal energies better for metallic systems than insulating 
ones. QMC, alternatively, provides an accurate descrip- 
tion of electron correlation^^ and can accurately de- 
scribe both states. We look for clues to the origin of the 
IMT in the DFT and QMC defect formation energies. 

In Fig. [3j the formation energy AEf of a Se defect is 
shown as a function of defect concentration. An indi- 
cation of a phase transition will appear as a kink (dis- 
continuity in the first or higher derivative) in the curve 
showing the total energy per atom vs. dopant concentra- 
tion; or equivalently as a kink in the curve showing defect 
formation energy vs. dopant concentration (since one is 
an affine transformation of the other). All experimentally 
known transitions in semiconductors currently appear to 
be continuous (notably a first-order Mott transition in 
LiCo02 is believed to exisi!^. In DFT, the defect forma- 
tion energy peaks in the vicinity of the IMT, although the 
discrete sampling renders it difficult to determine the or- 
der of any potential discontinuity On the insulating side, 
the defect formation energy increases with defect concen- 
tration, likely because correlation energy (a stabilizing 
contribution) tends to decrease with increasing electron 
density Thus the penalty for assigning additional Se 
atoms becomes more costly and ultimately renders the in- 
sulating state unstable. At higher Se concentrations than 
the IMT, the defect formation energy is relatively sta- 
ble, slightly decreasing with increasing defect concentra- 
tion. Using QMC (which, due to its computational cost 
was only performed for the Sei:Sii5, SeiiSiss, Sei:Sii27, 
and Sei:Si249, samples), the peak in AEf is amplified 
and shifted towards Sei:Sii27. This result remains in 
agreement with experiment, although QMC calculations 
should more accurately predict the transition point than 
DFT, which tends to overly delocalize electronic states. 
Also, the QMC correction to the defect formation energy 
is larger in the vicinity of the phase transition than for 
the metallic systems, indicating that as expected DFT 
more accurately simulates the metal than the insulator. 
This result suggests that many body effects are particu- 
larly important near the transition point and — although 
not a rigorous proof that the phase transition is driven 
by correlation — indicates that the IMT in Se:Si exhibits 
a strong Mott-like component. The differences between 
DFT and QMC descriptions of the Se:Si system clearly 
illustrate the importance of accurately treating electron 
correlation in the fundamental study of phase transitions. 

Finally, we comment on the implications of our 
analysis on the technological applications of silicon 
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hyperdoped with chalcogens, which has been considered 
as a candidate material for infrared absorbers and 
intermediate band photovoltaics (IBPVs). The anoma- 
lous sub-band gap absorption observed in hyperdoped 
samples can be understood in terms of the electronic 
band structure diagrams (showing a highly dispersive 
defect band that merges with the conduction bands) 
and the calculated absorption spectra. We believe that 
Se:Si is an intriguing candidate for an infrared detec- 
tor, exhibiting strong absorption down to low photon 
energies. Although the experimentally realized system 
may differ in important ways, the system we model 
computationally would perform poorly as an IBPV due 
to the lack of unoccupied states in the intermediate 
band, in agreement with Ref. 26, We note, however, 
that we have considered an uncompensated system 
in the above discussion, and that compensation may 
permit additional degrees of control over the properties 
of hyperdoped materials. We are optimistic, though. 



that we can use the framework described in this Letter 
to predict and screen new candidates for intermediate 
band properties. 
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